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ABSTRACT 

In this paper the problem of a finite plate containing col linear 
surface cracks is considered. The problem Is solved by using the line 
spring model with plane elasticity and Reissner's plate theory. The 
main purpose of the study is to investigate the effect of interaction 
between two cracks or between cracks and stress-free plate boundaries 
on the stress intensity factors and to provide extensive numerical 
results which may be useful in applications. First, some sample results 
are obtained and are compared with the existing finite element results,. 
Then the problem is solved for a single (internal) crack, two coll inear 
cracks and two corner cracks for wide range of relative dimensions. 
Particularly in corner cracks the agreement with the finite element 
solution is surprisingly very good. The results are obtained for semi- 
elliptic and rectangular crack profiles which may, in practice, corre- 
spond to two limiting cases of the actual profile of a subcrit ical ly 
growing surface crack. 

1 . Introduct ion 

Surface cracks are among the most common flaws in structural com- 
ponents, particularly in welded structures. Under cyclic loading or 
under static loading in the presence of corrosive environment any sur- 
face flaw has the potential of subcritical ly growing into a surface 
crack. Analysis of the structure containing such flaws is needed for 
modeling and prediction of the corresponding crack propagation rate. A 
review of the subject and a number of articles dealing with the analysis 
of the surface crack problem in plates may be found in [1]. At this 
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point the analytical treatment of the problem appears to be intractable. 
Therefore, the reliable solutions of the problem seem to be based oh 
numerical techniques, most notably on the finite element method (see, 
for example, [2] for the solution of a wide plate containing a semi- 
el liptic surface crack). In recent years, however, there has been some 
renewed interest in the application of the line spring model which 
was first described in [3] to the analysis of surface crack problems. 

The method was used in [4] in conjunction with Reissner's plate theory 
and the stress intensity factors for a semi-elliptic and a rectangular 
surface crack were calculated for a Wide plate under tension or bending. 
The semi-elliptic crack results described in [4] compare very favorably 
with the finite element solution given in [2]. 

In this paper the general problem is considered for a plate having 
a finite width. Analytically, it is known that if the stress fields 
Of more than one crack Or that of a crack and a stress-free boundary 
of the plate interact, there would be some magnification in the Stress 
intensity factors. The problem may therefore be important in plate 
structures having more than orte initial surface flaw or having a flaw 
near or at the boundary. Extensive finite element results for a single 
central or corner surface crack in a plate of finite width are given 
in [5] and [6] . Empirical ly developed express ions for stress intensity 
factors based on the results given irt [5] are also described in [?]. 

The present study was undertaken partly to show that the line spring 
model may be used for cracks in finite plates , part i cularly for corner 
cracks just as effectively as the infinite plate and partly tb supple- 
ment the results given in [5] and [6] by, for examp 1 e , considering the 
cases of a rectangular crack prof? le and col linear surface cracks. 

2. The General Formulation of the Problem 

The problem under consideration Is described in Fig. 1. It is 
assumed that x^x^ and x^x^ planes are planes of symmetry With respect 
to loading and geomet ry and the length of the plate in x 2 direction is 
relatively long compared to the width 2b so that in formulating the 

5 • 
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perturbation problem one may assume the plate to be infinitely long. 

Even though the numerical results are given for uniform tension in x ^ 
direction and cylindrical bending in x^ plane applied to the plate 
away from the crack region, as will be seen from the formulation of the 
problem, there is no restriction on the external loads provided in 
the absence of any cracks the membrane and bending resultants in x^ 
plane can be obtained for the given plate geometry and the applied loads. 

The problem is formulated for the collinear cracks shown in Fig. 1. 
The single central crack and the edge or the corner cracks are then 
considered as the special cases. One of the advantages of the line 
spring model is that the crack profile (as described by the function 
L(xj) giving the crack depth) can be arbitrary. However, the actual 
crack morphology studies indicate that for a given length 2a and a 
depth L q the crack profile may be bounded by a semi -ell ipse and a rec- 
tangle. Hence, in this paper the calculated results will be given only 
for these two limiting crack shapes. 

Ordinarily, the problems of in-plane loading (as expressed as a 
generalized plane stress problem) and bending of a plate are uncoupled. 
Consequently, the corresponding through crack problems can be solved 
independently. For the plate geometry shown in Fig. 1 the plane elasti- 
city and plate bending solutions are given in [8] and [9] , respectively. 
In the case of surface cracks, because of the absence of symmetry in 
thickness direction, the membrane and bending problems are clearly 
coupled. As in [9] in this paper, too, a transverse shear theory is used 
to formulate the bending component of the problem. The particular 
theory used is that of Reissner's [10] which is a sixth order theory 
and accounts for all three boundary conditions on the crack surfaces 
separately. 

Referring to Appendix A for normalized quantities and, for example, 
to [11] for the general formulation, the basic equations of the plate 
problem may be expressed as follows: 

V 4 (|> = 0 , 

V 4 w - 0 , 
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( 1 ) 

( 2 ) 



<V 2 i|) - - w = 0 , 

(3) 

k v 2 ft - n » 0 , 

(h) 

“xx ■ w (h,|,) • 0 yy * S* (h * ) 1 °xy * ' K3y <h * ) 

(5) 

**’%*¥*% * ’ 

(6) 

H xx ' iuT t 3 # + + 2 (, ‘ v> 3x3y 3 1 

(7) 

M yy = hP r V W" 7 < '" V) 3x3y 3 ’ 

(8) 

a(l-v) „ 3 2 * , x ,, , ,3-C _ 3 2 ni. 

M xy " IhT 1- P 3x3y + 2 0 v ’ { W ’ 

(9) 

V x " 3^ + 2 { \ v) 3y 3x ’ 

(10) 

. 3w . k n \ 3ft. . 3ij> 

V y = 3y '? (, - v) 3x + 3y 

(11) 


where, in the usual notation, F (or <|>) is the Airy stress function. 


N. . , 
IJ 


M. . , and V., ( i ,1=1 ,2) are the membrane, beiiding* and trahsverse 
ij i 


Shear resultants, 3-j and 3^ are the components of the rotation vector , 
u. , u_ and u, are the components of the d i spl acement vectbir , a is a 
1 ength parameter representing the crack Size (a =a for 0<e<d^b and a -d 


for c=0, d<b, Fig. 1), E and v are the elastic constants, the constants 
k and X are defined in Appendix A, \p and Q are auxi 1 iary functions 
defined in [1 1 ] , and the dimensions h, a, b, c, and d are shown in Fig. 1 . 

Because of symmetry, it is sufficient to consider the problem for 
0<Xj<b, 0<x 2 <» only. Thus, the membrane and bending problems of the 
plate must be solved Under the following boundary and symmetry conditions 
stated in terms of the normalized quantities (Fig. 1 and Appendix A)*. 


u(0,y) * 0, N (0,y) * 0, 0<y<» , 

Ay 


N xx (b',y) - 0, N xy (b' iY ) * 0, 0<y<» , 


-k- 
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04 ) 


N xy (x,0) ■ 0, 0<x<b ' , 

N yy (x,0) » 1 [-cr^x) + af(x)] , c'<x<d' , 05a) 

V(x,0) = 0, 0<x<c', d'<x<b' ; (15b) 

6 (0,y) * 0, M (0,y) - 0, V ( 0,y) » 0, 0<y<~ , 06) 

i x x y * 

M xx ( b '» y ) ° 0, M xy (b',y) » 0, V x (b',y) » 0, 0<y<« , 07) 

M (x,0) * 0, V (x, 0) = 0, 0<x<b' , 08) 

xy y 

M yy (x»0) ■ [-^(x) “ ^ * c ' <x<d ' , 09a) 

$ y (x,0) = 0, 0<x<c', d'<x<b' . (19b) 


The conditions stated above refer to the perturbation problem in which 

the crack surface tractions are the only nonzero external loads. Conse- 

quently, in addition to (12)-(19) it is required that 

N yy(x> <x> ) = °» N X y( x > w ) “ 0» 0£x<b' , (20) 

M yy (x,~) » 0, M xy (x,«>) = 0, V y (x,«) = 0, 0<x<b 1 . (21) 


The input functions and m^ which appear in (15a) and (19a) are 
defined by 

a «>(x) " N“ 2 ( x r 0)/h, m^fxj - 6M“ 2 (x 1 ,0)/h 2 (22) 

where n!\ (x 1 ,x 2 ) and mTj (x 1 ,x 2 ) , (I, >1,2) are the membrane and moment 
resultants in the plate under the actual applied loads in the absence 
of any cracks. The functions a(x) and m(x) are unknown and are defined 

by uU ) . 6m(x. ) 6M(a*x. 

(23) 


N (x. ) M / * v \ 6M(x.) 6M(a >c x) 

, .w . V k*— 


where the membrane load N(x^) and the bending moment M(x^) represent 
the stress component a 22 (x^,0,Xj) in the net ligament c<x^<d, 

-F<x 3 4-l. 
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In the bending problem the solution of the differential equations 
(2)-(4) satisfying the symmetry conditions (16) and the regularity 
conditions (21) may be expressed as follows [9]: 


w 


(x,y) * ~ | (A^+yA 2 ) e cosax dot 
o 

> * 

+ 7 I ( C i cosh ^ x + C 2 xs fnhgx) cos3y d6 * (24) 

o 

00 00 

( x *y) = \ | 8^6 1 sinax do + ^ J B 2 sinh r ? x sinBy d$ , (25) 


0 

oo 


cosax da 


^( x #y) * I" | t“A^(2Ka-y)A 2 ]e’ -ay 

o 

00 

+ “ | [-(C 1 +2icBC2)GOsh8x-C 2 xsinhex]cos3y dp, 

O 

where A. (a) , Bj (a) end C j ($) # ( J , 2 ) are unknown functions and 

r i - [ “ 2 + 7riW f >*i‘ t62 + mriW J 1 * 


(27) 


By substituting from (24)-(26) into (7), (9)*(ll) and by using five 
homogeneous conditions (17) and (18) five of the six unknown functions 
may be eliminated. The mixed boundary condition (19) would then deter- 
mine the sixth. 

Similarly from the plane stress solution of the plate satisfying 
the conditions (12), (14) and (20) the stresses and the y-componeht 
of the displacement may be expressed as [8] 

oo 

N XX (*.Y) “ “ “ | h^ (a) (1-oy)e cosaX da 
» o 

" | j [h 2 (B)cosh$x+6xhj (B) s i nhBx]cos8y dB , (28) 

o 
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Nyy ( x, y) 58 “ f" j h 1 (a)(I+ay)e" ay cosctx da 

o 

00 

+ -■ | [(h 2 +2hj)cosh8x+gxhjSlnh8x]cos6y d'6, (29) 

o 

00 

N xy^ X,y ^ " " f { a Y h i sinax da 

o 

00 

+ ~ | [(h 2 +hj)slnh8x + 8xhjCosh0x]sin$y d$ , (30) 

o 

00 h 

v(x,y) - ^ | (—£■ + ay)e~ ay cosax da 

o 

2 r h 2 i+K 

+ — I [(-j- + ~ 2 ~ h^)cosh$x + xh^sinh8x]sin8y d0 . (31) 

o 

In this case the unknown functions h^, h 2 and h^ are determined from 
the remaining boundary conditions (13) and (15). 

3. The Integral Equations 

If we now replace the mixed boundary conditions (15) and (19) 
respectively by 

— v(x,0) . - g 1 (x) , 0<x<b , (32) 

5“ 3 y (x,0) * g 2 (x) , 0<x<b , (33) 

it is seen that by using (17)* ( 1 8) , (13), (32) and (33) all nine 
unknown functions A { , B f , C., (i»l,2) and h. , (j*1,2,3) which appear 
in the formulation of the problem given in the previous section may be 
expressed in terms of the new unknown functions g^ and g 2 - From the 
definitions (32) and (33) It also follows that conditions (15b) and 
(19b) are equivalent to 
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( 34 ) 


g . (x) = 0 , (Kx<c* , d 1 <x<b ' , 0-1,2). , 
d' 

| g.(x)dx - 0 , (1*1,2) . (35) 


The functions g^ and g 2 may now be determined from the two remaining 
conditions (15a) and (19a). Referring to [8] and [9] for details, the 
following integral equations may be obtained from these two conditions: 

d ' 1 i cr^x) 

^F^x + t+x + k 1 (x,t) ‘ k^x.-tJlg^tJdt * — g~ , (36) 

c' 



m(x) 

~§r 


d' 

a*(i-v 2 ) { , r 3+v /_!_ + _U - 

2irhX 4 J 1 L 1+v 't~x t+x 1+ \> 

c' 



+ t^t c tir K 2<Y|t-x|) . ^ K 2 ( Y |t.x|)] + k 2 (x,t) 

m oo (x) 

- k 2 (x,-t)}g 2 (t)dt * - g g — , c ' <x<d ' , 


(37) 


where K 2 *is the modified Bessel function of the second kind, the 
Fredholm kernels k^x.t) and k 2 (x,t) are given in Appendix B and the 
constant y is given by 


h 

Y = 12(l-v ,4 )a* * 


(38) 


The functions a(x) and m(x) which appear in (36) and (37) are 


defined by (23) and represent the membrane and moment resultants of 
the tensile stress a 22 in the net 1 igament c'<x<d' . By using the plane 
strain solution for an edge crack occupying (h/2)-L<Xj<h/2 in a strip 
of thickness h (Fig. l) under membrane load N(xj) and bending, moment 


M(xj) (applied in x 2 x^ plane) and by expressing the rate of change of 
the potential energy in terms of crack closure energy and the change of 
compliance, a(x) and m(x) may be expressed in terms of the crack opening 
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displacement 2 v(x, 0 , 0 ) and the crack opening angle 2 $ (x, 0 ) as follows 
(see [l] and [4] for details): 


<t(x) - E [y tt (x)v(x) + Y tb (x)3 y (x)] , 
m(x) * 6 E [y bt (x)v(x) + Y bb (x)B (x)] , 


(39) 

(40) 


where the functions y.j, (i,j=t,b) depend on the local crack depth L(x) 
and hence are implicit functions of x. The algebraic expressions of 
these functions are given in [4] . From (32), ( 33 ) and (34) by observing 
that 


x x 

v(x,+ 0 ) - | g^ (t)dt , 3 y (x,+ 0 ) = f g 2 (t)dt , 


(41) 


and by using (39) and (40) , the integral equations ( 36 ) and ( 37 ) may 
then be expressed as 

x d' 

Y tt (x) j g 1 (t)dt "?Fj ^ th + tST + M** 1 * - k,(x.-t)] gi (t)dt 


X 

+ Y tb^{ 9 2 » c'<x<d' , 

c' 

Y bt (x){ 3 l (t)dt + y bb (x)J 9 2 <t)dt - l (-L. 


(42) 


+ \ \ _ ^ K (l~v) r 1 , 1 i , 4 r 1 v ( 1 . 1 \ 

t+x l+v ^ (t-x) d (t+x) 3 ■* + 1+v ^ t-x K 2 ^ t-x | ) 

+ ^(Ylt+x))] + k 2^ x,t: ^ " k 2^ x,_t ^ 9 2^ t ^ C * t 


meo(x) , c ' <x<d ' . 


(43) 


From the following asymptotic behavior of the Bessel function ^(z) for 
smal 1 values of z 
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K 2 (z) « pr - j + 0(z 2 1og z) (44) 

It can be shown that, as in (42), the Integral equation (43) has a 
simple Cauchy type singular kernel. We also note that the system of 
singular integral equations (42) and (43) must be solved under the 
additional conditions (35). 

After solving the integral equations (42) and (43) for g^ and g 2 
the Mode I stress intensity factor K at the leading edge of the crack 
may be obtained by substituting from (39)-(4l) into the following expres- 
sion giving K in a strip containing an edge crack of depth L and sub- 
jected to the membrane load a and bending moment m [4]: 

K(x) = 4T [cr(x)g t + m(x)g b ] (45) 

where g^. and g^ are functions of L/h and are obtained from the correspond- 
ing plane strain solution. From the results given in [12] the expressions 
for g t and g b valid in 0<L/h<0.8 may be obtained as follows: 


g t (s) = 47 (1.1216 + 6.5200s 2 - 12.3877s 4 + 89.0554s 6 

- I88.6080s 4 + 207.3870s 10 - 32.0524s 12 ) , (46a) 

g fa (s) = 47 (1.1202 - 1.8872s + 18. 0143s 2 - 87.3851s 3 

+ 241.9124s 4 - 319.9402s 5 + 168.0105s 6 ) , (46b) 


where s ■ L(x)/h. 

We now note that for 0<c'<d'<b the solution of the system of singu- 
lar integral equations is of the form 


9 j (*) 


G. (x) 

(x-c')*(d'-x)* 


, c'<x<d' 


0 - 1 , 2 ) , 


( 47 ) 


where the bounded unknown functions Gj and G 2 may easily be obtained 
by using the technique described, for example, in [ 13 ] . 
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The general crack geometry shown in Fig. 1 has two special cases. 
The first is the case of a symmetrically located single crack along 
-d'<x<d', (i.e., c lss 0, d'<b'). In this problem by using the symmetry 
considerations and by observing that g.(t) ■ -g . (**t ) , (1=1,2) , the 
integral equations (42) and (A 3 ) may be somewhat simplified as follows: 
x d' 

Y tt ( x > / 9 i (t > dt " 27 I C k, (x,t)]g 1 


-d 1 


-d 1 


X 

+ Y tb^j 9 2 ( t ) dt = jf a a> M » -d ,< x<d' » 

-d ' 

x x d' 

Y bt (x){ 9,(t)dt + Y bb (x)J g 2 (t)dt - [|^-p x 

-d' -d' -d' 


m 


TCT'^r K 2 <r|t ‘ x|) + k 2 (x,t)]g 2 ( t) d t 


1 

6E 


m (x) 
00 


-d ' <x<d 


(b 9 ) 


By using (44) it may again be shown that (49) has a simple Cauchy kernel 
and the solution of the integral equations is of the following form: 

F . (x) 

g, (x) = - r , -d'<x<d' , (i-1,2) . (50) 

1 (d'2-x 2 )* 


The second special case is that of corner cracks for which 0<c'<d'=b' 
In this case it may be shown that as x and t approach the end point b' 
simultaneously, the kernels k^ and k^ in (42) and (43) become unbounded. 

As shown in [8] and [ 9 ] the singular part of these kernels may be separ- 
ated and may be shown to be 


k ls (x» t) = k 2s (x,t) 


1 6(b'-x) . 4(b'-x) 2 

2b'-x-t Ub'-’x-t) 2 (2b'-x-t) J ’ 


where 


(51) 
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k } (x,t) » kj s (x * t ) + kj f (x,t) , (!»1 ,2) (52) 

s 

and k^ and k^ are bounded. Together with the Cauchy kernel 1/(t-x) in 

(42) and (43) • (51) constitutes a generalized Cauchy kernel* It may 

be observed that the generalized Cauchy kernel k (x,t) ■ 1/(t-x)+k. (x,t) 

9 ^5 

has the property that k g (x,b') = 0, k g ( b ' ^ t ) » 0 and consequently (t) 
and g 2 (t ) are nonsingular at t=b' [8] . Also, in this case the slrtgle- 
valuedness conditions (35) are not valid and, as pointed out in [8] , 
are not needed for a unique solution of the integral equations. 

4. The Results 

First, some sample problems are solved in order to compare the 
results obtained from the line spring model in this paper with that 
obtained from the finite element solutions given in [5] and [6] ; In [5] 
the single symmetric semi-elliptic surface crack problem is Considered 
for a finite plate under Uniform tension or cylindrical bending (i .e. , 
c=0, d<b, Fig. l). It is assumed that the half length of the plate is 
£=5d. Figures 2 and 3 show the cOmparisoh of the normalized stress 
intensity factors calculated along the crack front by the two methods. 

The normalizing stress intensity factor « N shown in these figures is 
defined by 

! K N = » k 38 ‘'l-^/d 2 (53) 

ahd is the stress intensity factor at the location x^ = 0* x^ * 0* x^ - 
L q , (lie., the end points of the minor axis) Of a flat elliptic craCk 
(with semi axes d and L q ) in an infinite solid subjected to Uniform 
tension cr 22 ■ cr^ In x 2 direction (c=0* Fig. 1). Note that, considering 
the simplicity of the line spring model, the agreement is not bad. One 

may also note that at the intersection point of the crack and the plate 
surface x » Xj/d ■ 1 the results based on the line spring model wOuld 
not be expected to be very good. Furthermore* at the singular point oh 
the free surface the power Of the stress singularity seems to be less 
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than 1/2 [14] .Hence, theoretically the stress intensity factor defined 
on the basis of conventional 1/2 power should tend to zero as the point 
on the crack front approaches the free surface at an angle of ir/2. Thus, 
strictly speaking, the bounded nonzero stress intensity factor given by 
the finite element solution at the surface do not seem to be correct either. 

Figures 4 and 5 show the comparison of the stress intensity fac- 
tors for a corner crack having the profile of a quarter ellipse and 
Obtained from the line spring model and the finite element solution given 
in [6] . It should be noted that the finite element results are obtained 
for a finite plate in which the half length is equal to the total width 
of the plate and the crack is only on one corner (see the insert in 
Fig. 4). However, since the crack length-to plate width ratio in both 
cases is relatively small (2a/2b = 1/10 in line spring and 2a/b = 1/5 
in finite element solution), the stress intensity factors for the two 
geometries should be approximately equal. The figures again show that 
the agreement is quite good. 

The calculated stress intensity factors are given in Tables 1-11. 

All stress intensity factors were calculated as a function of X s * x^/a*, 
(a*=d for a single crack, a*=a for two cracks. Fig. 1) defining the 

location along the crack front and of the relative dimensions of the 

< * 

Crack and the plate. The following notation and normal izing stress 
intensity factors are used in presenting the results: 

K 00 

CT b22^ r ’ 0 ’V “ * x = x i/ a * * (54) 

/2irr 
K (x) 

o.oo(r,0,x.) ■ - , x - x./a* (55) 

1 /2*r 1 

where supscripts b and t correspond to plates under bending and tension, 
respectively, 022 * s c 1 eava 9 e stress around the crack front, r and 0 
are the usual polar coordinates at the crack front in X 2 X^ plane (Fig. 1) 
and and are the corresponding Mode I stress intensity factors. 

The results are given for uniform membrane load l^ = N w and cylindrical 
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bending moment M 22 «* away from the crack region. The normalized 
stress intensity factors shown In the tables are defined by 


. Ku(x) 

K„(x) 


k b " ~K ^ — * k t " 

t 

K to ' 

(56) 

N 



K to ■ ( T> ^ 9 t <’o> > 

s o ’ L o /h • 

(57) 

6m 

ho a ^ g b (s o ) 

* s ■ L /h 
o o 

(58) 


where L q is the maximum crack depth and the functions g t and g b are 

given by (45) and (46). One may note that g. (s ) and g. (s ) are the 

T* O o O 

shape factors obtained from the corresponding plane strain solution of 
a plate with an edge crack of depth |» and, for the values of L /h shown 
in the tables, are given by [12] * 



Table 1 shows the normalized stress intensity factors at the deep- 
est penetration point of a centrally located single semielliptic surface 

crack (i.e., c*0, d<b. Fig, l) in a plate under uniform tension N or 

00 

bending M^. Here the crack profile is given by 



or 

L(x) - L o >W , (x * Xj/a* , a* = d) (60) 
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and hence x=0 is the deepest point on the crack front. This is also 
the point where k t assumes its maximum value. For b/h = 10 relatively 
complete and for other plate dimensions some sample results showing 
the variation of the stress intensity factors along the crack front 
are shown in Tables 2 and 3. Similar results are shown in Tables 4 
and 5 for a single surface crack with a rectangular profile (i.e., for 
L(x) = L q , - 1 <x<l ) . One may observe that, as expected, generally the 
stress intensity factors for the rectangular crack are higher than that 
for the semi-el 1 ipt ic crack. 

The results for two collinear semi -el 1 ipt ic surface cracks (Fig. 1) 
are shown in tables 6 and 7. Here the crack profile is defined by 
(Fig. 1) 

/ _ _ x r (c+a) 

L(x) » L o M-x 2 , x - — , -l<x<1 . (61) 


Table 6 shows the value k. (x*) , (i=b,t) and the location x - x* of the 
maximum stress intensity factor for various crack geometries in a plate 
for which b — 1 Oh and a = h. The factor 0 — a/(a+c) determines the crack 
location. Table 7 shows some sample results giving the distribution of 
the stress intensity factors along the crack front for two extreme crack 
locations considered. The skewness in this distribution does not seem 
to be very significant. 

The results for a plate containing two corner cracks having a pro- 
file of a quarter ellipse are shown in Tables 8 and 9 (Fig. 1). In 
this case the crack profile (or the crack depth) L is defined by 


L(x) = L / l-(-^) 2 
0 2 


x.-(c+a) 

x a , -1<X<1 

3 


( 62 ) 

i 


Table 8 shows the normalized Mode I stress intensity factors at the 
maximum penetration point of the crack which Is on the plate boundary 
x 58 b 1 (i.e., for x^ * b or x ** 1 or L * L q ) . Some results showing the 
distribution of the stress intensity factors are given in Table 9. The 
results were similar for al 1 crack geometries in that for plates under 
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tension and for those having shallow cracks under bending the maximum 
stress intensity factor was on the boundary x - b ', whereas for deep 
cracks in plates under bending K was maximum at the surface x-j = c or 
x = c' (Fig. 1). For corner cracks with a rectangular profile results 
similar to those shown in Tables 8 and 9 are given in Tables 10 and 11. 
For this crack geometry too one may note that generally the stress 
intensity factors for rectangular cracks are higher than those for the 
elliptic cracks. 

From the formulation of the problem it may be seen that all result? 
in the surface crack problem are dependent op the Poisson's ratio v 
of the plate. The stress intensity factors given in this paper are cal- 
culated for v » 0.3, However, as shown [9] » since the stress intensity 
factors are not very sensitive to the Poisson's ratio, the results 
given in Tables 1-11 should he valid for pearly all st ructural materials. 
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Appendix A 


The definition of normalized quantities 

x = x^/a* , y » x 2 /a* , z * x^a* , ( A .|) 

u - Uj/a* , v = u 2 /a* , w - u 3 /a* , ( A .2) 

♦ “ i^hE * e x = 0 i * 6 y - 3 2 ’ (A, 3) 

°XX = CT 11^ E » a yy = a 22 /E * ff xy * °12^ E » (A, k) 

N l . M.. 

N ctS = Tr * M a3 " iv^T » (a ’ 3) " ^X*y) » 0,j) * 0*2) , (A. 5) 

V x - V,/hB , V y = V 2 /hB , (a.6) 

B = f 2 0+v) .» K ^ I2(l-v 2 )a* 2 /h 2 . (A. 7) 

b' = b/a* , c' = c/a* , d' = d/a* 


In the problem described by Fig. 1, a* = a = (d-c)/2 for 0<c<d<b and 
a* ■ d for c = 0, d<b. 


- 18 - 



Appendix B 


The Fredholm kernels k 1 and k 2 which appear in the integral 
equations (36) and (37) 


kj (x,t) = | 


l+l» 3 b'e 


e " (2b '" t)6 r - 28 h' 

7 ^ 23 b'' — =W (-Cl + ( 3 + 23 b')e 6b ]coshpx 

* 6 — P 


- 2 Bxe 2Bb sinh 3 x-[ 23 x s inh 3 x+( 3 - 23 b ' 

+e 2Bb )cosh 3 x][l- 23 (b'-t)]}d 8 , (B.l) 

k 2 (x, t ) s(b'-t)] J±£^ re -' 21 > , - t -)e 

4 i -e 

o 

- I±£ 2r2X _ Br -(b'-t)B, -(b'-x)r 2 

1+v ^-2^' ^ B r 2 e ]e 


+ [(iL - Zb'B 2 H-e~ 2b,B ) ( - 2 g X) ±_ { 3 - 23 X) 

l M-v 1+v : -2b '8 M 1 1+v lKB U+e ’ 

I -e 


^xO-e-ZSx, . - 2 _ 6 (l+ e " 26x )}] 1 D] .-<a> ’-t-xlB 


+ 0je -(b'-x)6 e -(b'-t)r 2] . JK-B^O+e* 21 '**) ^. [0ie -(b'-t)B 


D Ll T 


* V- (b '- t)r 2 ] '-C.l e "(b'-x)r 2 }jg f 


1 -e 


2 b 'r 2 


(B. 2 ) 


23 _ /, _-2b'3\ l+e" 2b ' r 2 ^ . -2b 

y* r 2 (1 6 } 1 .;- 2 b T ~ r ' 2 ~ ° } 

. i+e" 2b, e P , 


[l-(b'-t)3]-(l-v)[| (b ' -t ) -<3 2 ]( 1 -e” 2b ' B ) , 


KY 


(B.3) 
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“2 - ■ 4 - <'- e ' 2b " S) ■ ^0-v)C-e- 2b '«, , 

1 “6 


D = Ab'B 2 e“ 2b ' 3 - + 2K8 3 )(l-e“ i,b ' 6 ) 


+ uZ « r 2 (1-e" 2b ' e ) 2 . 
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Table 1. The normalized stress intensity factors at the maximum 
penetration point (x=0) of a symmetrically located 
single semi-el 1 iptic surface crack in a plate under uni- 
form tension or bending (v=0.3). 


b 

■ 

L = 
o 

0 . 2h 

L = 
o 

0 . 4h 

L o ‘ 

0 . 6h 

L o ■ 

0 . 8h 

h 

■ 

k b ( o ) 

k t ( 0 ) 

k b ( o ) 

k t ( 0 ) 

k b (°) 

k t ( 0 ) 

k b ( o ) 

k t (0) 


0.5 

.709 

.729 

.308 

.390 

.0518 

.175 

-.0290 

0.0503 


0.6 

.737 

.755 

.342 

.421 

.0705 

.192 

-.0257 

.0555 


0.8 

• 777 

• 792 

.398 

.470 

.104 

.221 

-.0188 

.0648 


1 

.805 

.818 

.443 

.508 

.132 

.246 

-.0121 

.0730 

10 

4/3 

2 

.837 

.876 

.848 

.884 

.501 

.584 

.559 

.630 

.174 

.282 

-.0014 

.0848 

4 

.930 

.934 

.723 

.752 

.390 

.464 

.0726 

.155 


6 

• 953 

• 956 

.800 

.819 

.499 

.556 

.127 

.203 


8 

.967 

.969 

.853 

.865 

.592 

.634 

.190 

.256 


9-5 

.975 

.976 

.885 

.893 

.659 

.689 

.249 

.305 


9.61 

.976 

• 977 

.887 

.894 

.664 

.693 

.254 

.310 


9.8 

.977 

.978 

.891 

.898 

.672 

.700 

.264 

. 318 . 


0.5 

.709 

.729 


.390 

.0519 


-.0290 

.0503 


0.6 

.738 

.755 


.421 

.0706 


-.0256 

.0556 


0.8 

.778 

.792 

• 399 

.470 

.104 


-.0188 

.0649 

8 

1 

.805 . 

.818 

.444 

.509 

.133 

.247 

-.0120 

.0731 

2 

.877 

.885 

.586 

.632 

.246 

.341 

.0189 

.105 


4 

• 932 

.936 

.730 

.758 

.400 

.472 

.0774 

.159 


6 

.957 

.959 

.814 

.830 

.525 

.576 

.144 

.216 


7.69 

.971 

.972 

.867 

.876 

.626 

.660 

.223 

.282 


7.84 

.972 

.973 

.872 

• 880 

.635 

.667 

.233 

.290 


0.5 

.710 

.729 

L - " 

• 307 

.391 

.0521 


-.0289 

.0503 


0.6 

.738 

• 756 

.343 

.422 

.0710 

.192 

-.0256 

.0556 


0.9 

• 79 ^ 

.807 

.424 

.492 

j .122 

.235 

-.0152 

.0693 


1.2 

.827 

.839 

.483 

l .543 

.160 

.270 

-.0051 

.0807 

6 

1.5 

.851 

.861 

.530 

, .583 

i .196 

.301 

.0046 

.0910 


3 

.915 

.920 

.681 

• 715 

1 .341 

.423 

.0531 

.137 


4 

.930 

.934 

.723 

.752 

.390 

.464 

.0726 

.155 


5 

.953 

.955 

.802 

.818 

.507 

.560 

.136 

.208 


5.77 

.963 

.964 

.839 

.850 

.576 

.616 

.187 

.250 


5.88 

.964 

.965 

.844 

.855 

.587 

.625 

.197 

.258 
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Table 1 ( cont ) 


b 

d 

L = 
o 

0 . 2 h 

L .» 
o 

0 , 4 h 

L o = 

0 . 6 h 

L o = 

0 . 8 h 

h 

Hi 

k b (°) 

k t (o) 

k b (°) 

k t CO) 

k„Co) 

k t ( 0 ) 

k 6 t°) 

k t ( 0 ) 


0.5 

.711 

. 730 

.309 

.392 

.0528 

.176 

-.0289 

.0504 


0.666 

. 755 

.771 

.366 

.441 

.0839 

.204 

-.0231 

.0591 


0.8 

.780 

.795 

.403 

.474 

.106 

.223 

-.0184 

.0653 

2 i 

1 

.809 

.821 

.450 

.514 

• 137 

.250 

-.0112 

.0738 


1.33 

.843 

.853 

.512 

.568 

.183 

.289 

.0006 

.0866 


1.5 

.856 

.865 

.540 . 

.591 

.204 

.307 

.0068 

.0929 


2 

.886 

.893 

. 608 

.650 

.265 

.358 

.0257 

.111 


3.92 

.951 

.953 

.800 

.815 

.519 

.565 

.152 

.218 


0.5 

.716 

.735 


.398 

.0557 

.179 

-.0287 

.0508 


0.6 

.747 

.763 


.431 

.0768 

.197 

-.0249 

.0564 


0.8 

.791 

.804 

.421 

.488 

.117 

.232 

-.0166 

.0671 

2 

0.9 

.808 

.820 

.450 

.513 

.136 

.248 

-.0121 

.0122 

im 

1.0 

.823 

.843 

.477 

.537 

.156 

t .265 

-.0072 

. 071 b 


4/3 

.864 

.872 

.561 

.608 

, .224 

.321 

.0118 

.0961 


‘ 1.9 

.916 

.919 

.701 

.726 

• 385 

.450 

.0754 

.150 


1.96 

.920 

.924 

.718 


.411 

| 

.0903 

.162 
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Table 2. Distribution of the stress intensity factors along 
the crack front in a plate containing a single sym- 
metric semi-elliptic surface crack (b/h » 10, v * 
0.3, x'« Xj/d) . 



k b 

k t 

k b 

k t 

k b 

k t 

k b 

k t 

L o /h 

0 

.2 

0 

.4 

0 

.6 

0 

.8 

X 


b/h 

= 10 

d/h = 0.5 , 

v = 0 

.3 


0.929 

0.628 

.547 

.428 

.340 

.191 

HU 

.0486 

.444 

0.828 

.672 

.609 

.392 

.349 

.154 

HU 

.0314 

.472 

0 . 688 

.694 

.656 

.361 

.364 

.123 

.162 

.0113 

.510 

0.516 

.704 

.691 

• 336 

.376 


.169 

-.0061 

.512 

0.319 

.708 

.715 

.318 

.385 


.173 

-.0187 


0.108 

.709 

.727 

.308 

.390 

.0535 

.175 

-.0276 

mm 

0 

• 709 

.729 

• 307 

.390 

.0518 

.175 

-.0290 

.503 


b/h = 10 , d/h = 1 , v = 0.3 


0.929 

.631 

.545 

.505 

• 391 

.272 

.205 

.0809 

.0649 

0.828 

.709 

.639 

.496 

.426 

.239 

.215 

.0621 

.0677 

0.688 

.756 

.710 

.480 

.457 

.209 

.226 

.0396 

.0718 

0.516 

.783 

.762 

.464 

.482 

.177 

.236 

.0183 

.0729 

0.319 

.798 

.798 

.451 

• 499 

.149 

.243 

.0163 

.0724 

0.108 

.804 

.816 

.444 

.507 

.134 

.246 

-.0103 

.0728 

0 

.805 

.818 

.443 

.508 

.132 

1 .246 

-.0121 

j .0730 


b/h = 10 , d/h = 4 , v = 0.3 


0.929 

.623 

.535 

.561 

.420 

.402 

.285 

.168 

.121 

0.828 

.739 

.661 

.626 

.517 

.420 

.339 

.163 

.137 

0.688 

.819 

• 763 

.666 

.601 

.426 

• 387 

.144 

.150 

0.516 

.875 

.844 

.695 

.671 

.418 

.425 

.120 

. 1 56 

0.319 

.910 

.901 

.713 

.722 

.402 

.451 

.0953 

.156 

0.108 

.927 

• 930 

.722 

.748 

.391 

.463 

.0756 

.155 

0 

.930 

.934 

.723 

.752 

.390 

.464 

.0726 

.155 


b/h = 10 , d/h = 8 , v = 0.3 


0.929 

.622 

.533 

.571 

.423 

.453 

.316 

.238 

.170 

0.828 

.747 

.667 

.665 

.542 

.513 

.403 

.260 

.209 

0.688 

.837 

.778 

.735 

.653 

.560 

.487 

.261 

.240 

0.516 

.901 

.868 

.791 

.749 

.586 

.558 

.245 

.256 

0.319 

.944 

.931 

.830 

.821 

.593 

.607 

.219 

.259 

0.108 

.965 

.965 

.850 

.860 

.592 

i .631 

.194 

.256 

0 

.967 

.969 

.853 

.865 

.592 

.634 

.190 

.256 


































Table 2 (cont.) 


I 

» 

I 

t 

r 

I 

1 



mm 

mm 

is 

is 

mm 

IS 

k b 

k t 

L o /h 

0 , 

.2 

0. 

4 

0.6 

0. 

.8 

X 


b/h 

- 10 

d/h * 

9.8 , 

v = 0.3 


0.929 

.629 

.538 

.597 

.442 

.508 

.355 

.312 

.225 

0.828 

.753 

.673 


.562 

.572 

.446 

.341 

.270 

0.688 

.844 

.784 

.763 

.675 

.626 

.536 

.345 

.305 

0.516 

.909 

.875 


.775 

.658 

.614 

.328 

.323 

0.319 

.952 

.939 

.865 

.851 

.670 

.669 

.298 

.324 

0.108 

.974 

.973 

.888 

.892 

.672 

.697 

.268 

.319 

0 

.977 

.978 

.891 

.898 

.672 

.700 

.264 

.318 
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Table 3. Distribution of the stress intensity factors 

along the crack front in a plate containing a single 
symmetric semi-elliptic surface crack (b/h * 2,4.6: 
v=0.3). 



k b 

k t 

k b 

k t 

k b 

k , 

k b 

k t 

L /h 
o 

0 

.2 

0.4 

0 

.6 

0 

.8 

X 



b/h 

- 2, d/h = 

1, v = 

0.3 


0.929 

.646 

.559 

.542 

.421 

.306 

.232 

.0941 

0.0752 

0.828 

.726 

.654 

.533 

.456 

.271 

.240 

.0736 

.0768 

0.688 

.774 

.726 

.517 

.487 

.238 

.249 

.0492 

.0796 

0.516 

.801 

.779 

.500 

• 511 

.204 

.257 

.0261 

.0793 

0.319 

.816 

.814 

.486 

• 527 

.174 

.262 

.0077 

.0777 

0.108 

.823 

.832 

.478 

.536 

.158 

.264 

-.0053 

.0774 

0 

.823 

.834 

.477 

• 537 

.156 

.265 

-.0072 

.0774 


b/h = 4, d/h = 1 , v = 0.3 


0.929 

.634 

.548 

• 512 

.397 

.278 

.210 

.0833 

.0668 

0.828 

• 713 

.642 

.504 

.432 

.245 

.220 

.0642 

.0694 

0.688 

.760 

.713 

.488 

.463 

.214 

.230 

.0414 

.0733 

0.516 

.787 

.766 

.471 

.488 

.182 

.240 

.0198 

.0741 

0.319 

.802 

.801 

.458 

.505 

.154 

.246 

.0028 

.0734 

0.108 

.808 

.819 

.451 

.513 

.139 

.249 

-.0094 

.0737 

0 

.809 

.821 

.450 

.514 

.137 

.250 

-.0112 

.0738 


b/h = 6, d/h = 1.2, v - 0.3 


0.929 

.632 

.545 

.522 

.402 

.296 

.221 

.0921 

.0723 

0.828 

.717 

.645 

.523 

.446 

. 266 

.234 

.0732 

.0754 

0.688 

.770 

.722 

.513 

.483 

.237 

.247 

.0501 

.0796 

0.516 

.801 

• 778 

.501 

.512 

.206 

.258 

.0277 

.0808 

0.319 

.819 

.817 

.490 

.532 

.178 

. 266 

.0096 

.0802 

0.108 

.827 

.836 

.474 

.541 

.162 

.270 

-.0032 

.0806 

0 

.827 

.839 

.483 

.543 

.160 

.270 

-.0051 

.0807 



Table 4. The normalized stress intensity factors at the center 
(x®0) of a single symmetric rectangular surface crack 
in a plate under tension or bending (v = 0.3) . 


b 

d 

L o “ 

0.2h 

L o- 

0.4h 

<0“ 

0.6h 

L : = 
0 

0.8h 

h 


k b (o) 

k t (0) 

k b (Q) 

k t (o) 

k b (0) 

k t (0) 

yo) 

k t (Q) 


0.5 

.765 

.784 

.340 

.429 

.0607 

.194 

-.0316 

.0599 

10 

2 

• 915 

.922 

.652 

.699 

.284 

.388 

.0261 

.122 

5 

• 970 

.973 

.847 

.868 

.544 

.611 

• 134 

.222 



.999 

.999 

.987 

.989 

• 914 

• 927 

.557 

• 603 


B9I 

.766 

• 785 

.340 

.429 

.0608 

.194 

wan 

.0599 

8 

mm 

.853 

.865 

.496 

• 563 

.154 

.276 

Krctfl 

.0851 

1 1 

.963 

.966 

.814 

.840 

.487 

.562 

.104 

.195 


Bt9 

• 998 

• 998 

.982 

• 985 

.892 

.907 

.503 

.554 


0.5 

.766 

.785 

.341 

.429 

0.0610 

.194 

-.031$ 

.Q$00 

6 

1 

.855 

.867 

.498 

.566 

.155 

.277 

-.0103 

.0854 


3 

.951 

• 955 

.767 

.797 

.414 

.500 

.0721 

.165 


5.88 

.997 

• 998 

.975 

.978 

.857 

.878 

.434 

.491 


0.5 

.768 

.787 

.343 

• 431 

0.0619 

.195 

-.0315 


4 

1 

.859 

.870 

• 505 

• 571 

.159 

.281 

-.0095 

.0863 


2 

.930 

.936 

.690 

.732 

.320 

.419 

.0370 

• 133 


3.92 

.996 

.996 

• 959 

.965 

• 797 

.826 

•341 

.408 



.776 

.794 

.352 

.439 

.0655 

.198 

-.0312 

.0609 

2 

Ik 

.880 

.890 

.545 

.606 

.186 

.304 

-.0041 

.0923 



.9^1 

• 945 

.710 

.749 

• 334 

.432 

.0395 

.135 

a 

19 

• 990 

.991 

.916 

.927 

.666. 

.715 

.205 

.285 















































Table 5. Distribution of the stress intensity factors along 
the crack front in a plate containing a single sym- 
metric rectangular surface crack, x = x^/d. 



Hi 

IH 

■■ 

a 

HI 

Hi 

I 

k t 

V h 

0 

.2 

0 

.4 

0 

.6 

0 

.8 

X 



b/h = 

2, d/h = 1, V 

= 0.3 



0.929 

.585 

.618 

.233 

.334 


.159 

-.0295 

0.0458 

0.828 

.737 

.759 

.354 

.440 

.0798 

.209 

-.0261 

.0619 

0.688 

.814 

.829 

.439 

.514 

.122 

.248 

-.0190 

.0741 

0.516 

.852 

.864 

.495 

.562 

.154 

.276 

-.0120 

.0831 

0.319 

.871 

.881 

.528 

• 591 

.174 

.294 

-.0070 

.0890 

0.108 

.879 

.889 

.543 

.605 

.184 

.302 

-.0044 

.0920 

0 

.880 

.890 

.545 

.606 

.186 

.304 

-.0041 

.0923 




b/h = 

6, d/h = 1 , v 

= 0.3 



0.929 

.566 

.601 

.210 

.314 

.0181 

.149 

-.0302 

BiBIl 

0.828 

• 715 

.738 

.321 

.411 

.0623 

.194 

-.0283 


0.688 

.789 

.806 

• 399 

.480 

.0996 

.228 

-.0227 


0.516 


.841 

.451 

.524 

.127 

.253 

-.0169 

llltM 9 

0.319 

.846 

.858 

.482 

.551 

.145 

.269 

-.0127 

.0825 

0.108 

.854 

.866 

.496 

.564 

.154 

.276 

-.0105 

.0851 

0 

.855 

.867 

.498 

.566 

.155 

.277 

-.0103 

.0854 




b/h = 

10, d/h = 1 , V = 0.; 

5 


0.929 

.423 

.470 

.112 


-.0172 


! 

-.0293 

.0309 

0.828 

.574 

.609 

• 191 

.298 

.0038 

Bra 

-.0343 

.0417 

0.688 

.667 

.694 

.252 

.352 

.0250 

Bra 

-.0350 

.0492 

0.516 

.721 

.744 

.297 

.390 

.0421 

.177 

-.0339 

.0545 

0.319 

• 751 

.771 


.415 

.0539 

.188 

-.0325 

.0580 

0. 108 

.764 

.783 


.427 

.0599 

.193 

-.0317 

.0597 

0 

.765 

.784 


.429 

.0607 

.194 

-.0316 

.0599 
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Table 6. The location x=x* and magnitude k[j(x*) and k t (x*) of the 
normalized stress intensity factors in a plate containing 
two col linear semi -elliptic surface cracks, D=a/(a+c). 


D 


L o ' 

0.2h 

L - 
0 

0.4h 

L = 
o 

0.6h 

L - 
o 

0.8h 

k b ( **> 

k t (x*) 

k b (x*) 

k t (x*) 

k b (x*) 

k t (x*) 

k b (x*) 

k t (x*) 

0.112 

X* 

0.2 

0.05 

.929 

.319 

.929 

.929 

.929 

.929 


k(x*) 

.831 

• 839 

.649 

.554 

.409 

.308 

.138 

.107 

0.125 

X* 

0 

0 

.929 

.040 

• 929 

.108 

.929 

.516 


k(x*) 

.812 

.824 

.522 

.518 

.287 

• 523 

.867 

.756 

.250 

X* 

0 

0 

+.929 

0 

.929 

0 

• 929 

+ .516 


k (x*) 

.807 

.820 

.509 

.512 

.275 

.248 

.0822 

.735 

0.5 

x* 

0 

0 

-.929 

0 

-.929 

0 

-.828 

.516 


k(x*) 

.81 1 

.823 

.521 

• 517 

.285 

.251 

.0858 

.0744 

0.75 

x* 

-0.50 

0 

-.929 

-.050 

-.929 

-.108 

-.929 

-.688 


k(x*) 

.818 

.829 

.550 

.528 

.310 

.259 

.0951 

.786 
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Table 7» Distribution of the normalized stress intensity 

factors along the crack front in a plate contain- 
ing two coll inear semi-elliptic surface cracks, 
x = [x 1 ~(c+a)]/a (Fig. 1). 


* 


* 



k b 

k t 

k b 

k t 

k b 

k t 

k b 

k t 

L o /h 

0, 

.2 

0, 

.4 

0.6 

0, 

.8 

X 

b/h 

= 10, d/h = 1 , 

, D = a/(c+a) = 

0.112. 

, v=0.3 


0.929 

.688 

.596 

.649 

.505 

.409 

.308 

.138 

.107 

0.828 

.766 

.689 

.623 

.527 

.351 

.300 

.106 

.102 

0.688 

.805 

.754 

.584 

.541 

.297 

.294 

.0720 

.0975 

0.516 

.824 

.798 

.548 

.550 

.246 

.289 

.0411 

.0915 

0.319 

.831 

.827 

.519 

.554 

.204 

.285 

.0175 

.0858 

0.108 

.831 

.839 

.500 

.553 

.178 

.280 

.0013 

.0829 

0 

.829 

.839 

.494 

.550 

• 173 

.278 

-.0016 

.0821 

0.108 

.826 

.835 

.491 

.546 

.172 

.275 

-.0004 

.0814 

■0.319 

,.816 

.814 

.492 

.532 

.184 

.269 

.0117 

.0809 

•0.516 

.799 

• 776 

.500 

.512 

.209 

.261 

.0293 

.0818 

0.688 

1.769 

.721 

.513 

.484 

.240 

.250 

.0516 

.0814 

•0.828 

.720 

.649 


KBl 

.270 

.239 

.0751 

.0780 

•0.929 

.640 

.553 


Baa 

.303 

.229 

.0949 

.0758 


b/h =10, d/h = 1 , D = a/(c+a)=0. 75, v = 0.3 


0.929 

.637 

.551 

.521 

.404 

.288 

.217 

.0872 

.0698 

0.828 

.716 

.645 

.514 

.440 

.254 

.227 

.0678 

.0721 

0.688 

.764 

• 717 

.499 

.472 

.224 

.237 

.0446 

.0757 

0.516 

.793 

.771 

.484 

.498 

.192 

.247 

.0225 

.0763 

0.319 

.809 

.807 

.472 

.516 

.164 

.254 

.0050 

.0753 

0.108 

.816 

.826 

.467 

.526 

.149 

.258 

-.0075 

.0754 

0 

.818 


.467 

.528 

.148 

.258 

-.0093 

.0755 

•0.108 

.818 


.469 


.151 

.259 

-.0073 

.0755 

■0.319 

.814 

.812 

.480 


.169 

.258 

.0057 

.0760 

•0.516 

.801 

.778 

.497 

.509 

.200 

.253 

.0243 

.0778 

•0.688 


• 727 

.517 

.488 

.236 

.247 

.0481 

.0786 

•0.828 


.657 

.538 

.460 

.272 

.240 

.0735 

.0766 

•0.929 


.563 

.550 

.427 

.310 

.235 

.0951 

.0757 . 
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Table 8. The normalized intensity factors on the edges (x^+b') of 
a plate containing two symmetric corner cracks having a 
prof i le of a quarter el 1 ipse (Fig. 1 ) . 


b 

n 

mm 



"Oh 

L o " 

"Oh 

L = 
0 

0 . 8 h 

h 

u 

Si 

nan 

BS8I 

ngn 

naan 

|Dig« 

nan 

mam 


0.25 

.775 

• 790 

.380 

.485 

.0975 

.219 

-.0172 

.0678 


0.3 

.797 

.810 

.415 

.485 

.120 

.239 

-.0117 

.0743 


0.4 

.828 

.840 

.4 73 

.535 

.159 

.271 

-.0014 

.0857 

2 

eh 

.852 

.862 

.522 

.477 

.197 

.303 

.0089 

.0963 


EH 

.872 

.880 

.568 

.616 

.234 

.334 

.0199 

.107 


0.7 

.889 

.896 

.610 

.652 

.273 

.366 

.0321 

.118 


0.8 

.905 

.910 

.653 

.688 

.317 

.401 

.0470 

.131 


0.26 

.777 

.792 

.384 

.459 

.102 

.223 

-.0152 

.0700 


0 .4 

.821 

.833 

.463 

.527 

.156 

.269 

.0003 

.0873 


0.6 

.858 

.867 

.539 

.593 

.215 

.319 

.0179 

.106 

4 

0.8 

.883 

.890 

.597 

.642 

.264 

.361 

.0336 

.121 


1 

.901 

• 907 

.644 

.683 

.310 

.399 

.0492 

.136 


1 .2 

.916 

.921 

.685 

.718 

.354 

.435 

.0657 

.150 


1 .4 

.929 

• 933 

.722 

.750 

.398 

.471 

.0838 

.166 


1.6 

.939 

.942 

.756 

.779 

.443 

.508 

.105 

.184 


0.27 

.781 

• 796 

.391 

.464 

.106 

.226 

-.0140 

.0714 


0.3 

.792 

.806 

.410 

.481 

.119 

.237 

-.0105 

.0755 


0.6 

.856 

.866 

.536 

.591 

.214 

.319 

.0187 

.107 

6 

0.9 

.889 

.896 

.613 

.657 

t .281 

.376 

.0409 

.128 


1.2 

.910 

.916 

.669 

• 705 

.337 

.422 

.0609 

.147 


1.5 

.926 

.930 

.713 

.744 

1 .387 

.464 

.0809 

.165 


1.8 

.938 

.941 

.750 

.776 

.434 

j .503 

i .102 

.183 


2.1 

.948 

.950 

.782 

.803 

.479 

.540 

.124 

.202 


2.4 

• 956 

.958 

.811 

.828 

.523 

.576 

.149 

.223 


.28 

.785 

.799 

.397 

.470 

.110 

.230 

-.0129 

.0727 


0.4 

.821 

.833 

.462 

.526 

.156 

.269 

.0004 

.0875 


0.8 

.879 

.887 

.589 

.636 

.260 

.358 

.0343 

.122 

8 

1 .2 

.908 

.914 

. 665 

.702 

.334 

.420 

.0607 

.147 

1.6 

.927 

.932 

.718 

.748 

.394 

.470 

.0846 

.168 


2 

.941 

.944 

.760 

.784 

.447 

.514 

.108 

.189 


2.4 

.951 

.953 

.793 

.813 

.494 

.554 

.133 

.210 


2.8 

• 959 

.961 

.821 

.837 

.538 

.591 

.158 

.232 


3.2 

.965 

.967 

.845 

.858 

.580 

.626 

.186 

.255 
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Table 8 - cont 


b 

H 




EH 

L 0 . 6 h 

O 

L = 0 . 8 h 
o 

h 

m 



HS 

BSH 

Bsai 


BHI 

k t ( b ') 


0.25 

.772 

.787 

.376 

.452 

.0967 

.218 

-.0165 

.0684 


0.75 

.873 

.882 

.576 

.625 

.249 

.349 

.0307 

.118 


1 

.895 

.902 

.630 

.672 

.299 

.391 

.0483 

.135 


1.5 

.922 

• 927 

.704 

• 736 

.378 

.457 

.0786 

.163 

10 

2 

.939 

.943 

.755 

.780 

.440 

.510 

.106 

.188 

2.5 

.951 

.953 

.793 

.813 

.494 

.555 

.133 

.211 


3.0 

.959 

.961 

.823 

.840 

.541 

.594 

.160 

.234 


3-5 

.966 

.968 

.848 

.861 

.584 

.630 

.188 

.258 


4.0 

• 972 

.973 

.869 

.880 

.624 

.664 

.218 

.283 


4.5 

.976 

.977 

.888 

.896 

.664 

.696 

.252 

.312 

20 

m 

.895 

.901 

.629 

.671 

.298 

• 390 

.0483 

.135 
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Table 9. Distribution of the normalized stress intensity 

factors along the crack front [n a plate containing 
two (elliptic) corner cracks, x * [x.-(c+a)]/a 
(Fig. 1). 



E9 

■a 

k b 

m 

k b 

m 

k b 

k t 

L o /h 

0 

.2 

0.4 

0 

.6 

0 

.8 

_ 

X 



b/h = 

2, a/h = 0.5, v = 

0.3 


0.999 

.852 

.862 


.577 

.197 

.303 

.0089 

.0963 

0.936 

.846 

.856 

.515 

.571 

.191 

.297 

.0073 

.0936 

0.784 

.834 

.843 

.503 

.557 

.182 

.286 

.0050 

.0883 

0.558 

.824 

.828 

• 493 

.543 

.177 

.274 

.0064 

.0834 

0.279 

.813 

* 808 

.492 


.184 

.266 

.0138 

.0805 

-0.026 

.799 

.777 

.498 

.510 

.204 

.257 

.0263 

• 0798 

-0.329 

.776 

.732 

.511 

.488 

.231 

.248 

.0450 

.0794 

-0.600 

.736 

.669 

.526 

.460 

.261 

.240 

.0679 

.0768 

-0.815 

.668 

.583 

.537 

.427 

.294 

.231 

.0882 

.0748 

-0.953 

.549 

.460 

.532 

• 390 



• 336 

.232 

.112 

.0772 




b/h = 

8, a/h = 0.8, v = 

0.3 


0.999 

.879 

.887 

.589 

• 636 

.260 

• 358 

.0343 

.122 

0.936 

.874 

.882 

.582 

.630 

.253 

• 351 

.0317 

.118 

0.784 

.866 

.872 

.570 

.617 

.242 

• 339 

.0277 

.112 

0.558 

.857 

.859 

.561 

.602 

.237 


.0281 

.105 

0.279 

.844 

.836 

.557 

• 583 

.243 

.314 

.0357 

.101 

-0.026 

.825 

.800 

.557 

.559 

.259 

.302 

.0489 

.0990 

-0.329 

j .793 

.746 

.558 

• 525 

.282 

.286 

.0677 

.0966 

-0.600 

.741 

.671 

.558 

.481 

.304 

.269 

.0896 

.0917 

-0.815 

1.658 

• 573 

.547 

.428 

.326 

.249 

.108 

.0869 

-0.953 

; .521 

.434 

.502 

.361 

-345 

.231 

! . 127 

.0850 
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Table 10. The normalized stress intensity factors at the edges x = +b' 

of a plate containing two symmetric rectangular corner cracks. 


b 


L - i 
o 

0.2h 

L = 
0 

0.4h 

L o = 

0.6h 

L - 1 
0 

0.8h 

h 

m 

k b (b ') 

k t (b') 



k b (b ') 

k t (b') 

k b< b '> 

Mb') 


0.25 

.821 

.835 

.415 

.494 

.108 


-.0185 


2 

0.5 

.895 

.903 

.581 

.638 

.223 


-.0119' 

Bffii 


0.8 

.954 

.958 

.754 

.787 

.388 

.**77 

.0620 

■ is 


0.26 

.820 

.835 

.419 

.497 

.112 

.242 

-.0163 

.0797 

4 

0.4 

.860 

.871 

.507 

.574 

.174 

.295 

.0014 

.0985 


l 

.937 

.942 

.716 

.755 

• 359 

.**53 

.0595 

.154 


1 .6 

.976 

.978 

.856 

.876 

.550 

.617 

.139 

.221 


0.27 

.823 

.838 

.426 

wem 

.117 

.246 

-.0149 

0.0812 

6 

0.6 

.891 

.900 

.589 

KS-Jb 

! .240 

.353 

.0227 

.120 


1.6 

.956 

.960 

.788 

MEumI! 

.453 

.534 

.0983 

. 190 


2.4 

.984 

.985 

.902 

1 


.648 

.700 

.202 

.283 


0.28 

.827 

.841 

.433 

.510 

.122 

.250 

-.0136 

.0826 

8 

0.8 

.912 

.919 

.648 

.696 

.294 

• 399 

.0409 

.138 

2 

• 967 

.970 

.833 

.856 

.525 

.595 

.133 

.222 


3.2 

.988 

.989 

• 927 

.937 

.714 

.756 

.255 

.331 


0.27 

.823 

.837 

.425 

.503 

Bum 

.246 

-.0150 

.0811 

10 

1 

.927 

• 933 

.692 

.734 


.439 

.0573 

• 153 


2.5 

.974 

.976 

.864 

.882 


.643 

.165 

.251 


4 

• 991 

• 992 

.943 

.951 


. 796 

.302 

• 374 


- 33 - 












































Table 11. Distribution of the normalized stress intensity 

factors in a plate with rectangular corner cracks, 
x = kj-(c+a)j/a. 



is 

IB 

HU 

k t 

IB! 

IB 1 

■B! 

k t 

L o /h 

0 . 

2 

0 . 

4 

0 . 

6 

0 . 

8 

X 


CM 

11 

sz 

-Q 

a/h = 0 . 5 , 

v - 0 . 

3 


0.999 

.895 

.903 

.581 

.638 

.223 

• 337 

.0119 

.109 

0.936 

.892 

.901 

.576 

.634 

.218 

.332 

.0102 

.107 

0.784 

.887 

.896 

.564 

.623 

.207 

.323 

.0059 

.103 

0.558 

.879 

.889 

.548 

.609 

.193 

.310 

.0003 

.0967 

0.279 

.868 

.879 

.525 

.589 

.175 

.295 

-.0057 

.0902 

- 0.026 

.851 

.863 

.493 

.561 

.153 

.275 

-.0122 

.0828 

- 0.329 

.818 

.833 

.444 

.518 

.124 

.249 

-.0195 

.0739 

- 0.600 

.756 

.776 

.370 

.454 

.0852 

.214 

-.0268 

.0626 

- 0.815 

.630 

. 660 

.262 

.359 

.0373 

.168 

-.0314 

.0481 

- 0.953 

.385 

.434 

.115 

.229 

-.0098 

.107 

-.0260 

.0284 


b/h - 8 , a/h = 0 . 8 , v = 0.3 


0.999 

.912 

• 919 

.648 

.696 

.294 

.399 

.0409 

.138 

0.936 

.911 

.918 

.643 

.692 

.289 

.394 

.0383 

.135 

0.784 

.907 

.915 

.633 

.683 

.277 

.384 

.0322 

.129 

0.558 

! -902 

.910 

.619 

.671 

.261 

.370 

.0249 

.121 

0.279 

.894 

.902 

.598 

.652 

.242 

.353 

.0172 

.114 

- 0.026 

. 880 

.890 

.567 

.625 

.217 

• 331 

.0086 

.105 

- 0.329 

.856 

.868 

.520 

.584 

.183 

.301 

-.0015 

.0934 

- 0.600 

.809 

.824 

.446 

.520 

.137 

.260 

-.0133 

.0792 

- 0.815 

.705 

.729 

.331 

.420 

.0759 

.205 

-.0245 

.0610 

- 0.953 

.462 

.505 

.160 


.0073 

.131 

-.0263 

.0367 
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x * x | / d 


Comparison of stress intensity factors calculated by the 
finite element and line spring methods in a plate containing 
a symmetrically located semi -el 1 iptic surface crack and sub- 
jected to uniform tension. 





Line Spring 


x = x , /d 


0 . 2 ,- 


Lme Spring 
R 


x = x,/d 


Compar-Uon of the stress intensity factors calculated by the 
finite element and line spring methods in a plate containing 
a single syinnetric semi-elliptic surface crack and subjected 
to uniform bending. 






Lo 

2a 


0.4 



1 


2a 


T77 


L. / ^ 






Fig- 4 Comparison of the stress intensity factors calculated by the 
finite element and line spring methods in a plate containing 
elliptic corner cracks and subjected to uniform tension, 
L 0 /2a = 0.4, b/2a = 5. 
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